DAMTP-2012-49 
LMU-ASC 49/12 



On the Hamiltonian form of 3D massive gravity 



C/3 



Olaf Hohm 1 , Alasdair Routh 2 , Paul K. Townsend 2 
and Baocheng Zhang 2 ' 3 



1 Arnold- Sommerf eld- Center for Theoretical Physics, Fakultat fur Physik, 



Ludwig- Maximilians- Universitat Miinchen, Theresienstrasse 37, 80333 Miinchen, 



Germany 

email: olaf . hohm@physik . uni-muenchen . de 



2 Department of Applied Mathematics and Theoretical Physics, 
Centre for Mathematical Sciences, University of Cambridge, 



Wilberforce Road, Cambridge, CBS OWA, U.K. 
email: A.J.Routh@damtp.cam.ac.uk, P.K.Townsend@damtp.cam.ac.uk 

3 State Key Laboratory of Magnetic Resonances and Atomic and Molecular Physics, 
Wuhan Institute of Physics and Mathematics, 
Chinese Academy of Sciences, 
Wuhan 430071, China 

e-mail: zhangbc@wipm .ac.cn 
ABSTRACT 

We present a "Chern-Simons-like" action for the "general massive gravity" model 
propagating two spin-2 modes with independent masses in three spacetime dimensions 
(3D), and we use it to find a simple Hamiltonian form of this model. The number of local 
degrees of freedom, determined by the dimension of the physical phase space, agrees with 
a linearized analysis except in some limits, in particular that yielding "topologically new 
massive gravity" , which therefore suffers from a linearization instability 



1 Introduction 

Massive gravity models have been intensively investigated over the past few years, partly 
motivated by the idea that some infra-red modification of General Relativity (GR) could 
provide an alternative to dark energy; see e.g. [1] for an overview. As for many other 
issues in GR, it is useful to consider how things simplify in the context of a three- 
dimensional (3D) spacetime. Although the long-standing "topologically massive gravity" 
(TMG) model [2] shows that a massive graviton is possible in 3D, this is achieved by 
the introduction of the parity violating, Lorentz-Chern-Simons (LCS) term, which is 
the action for 3D conformal gravity [3j. In contrast, the more recent "new massive 
gravity" (NMG) model [HE] achieves a similar effect without breaking parity through 
the introduction of a particular curvature-squared term; this model exhibits in simplified 
form some of the features of ghost-free higher- dimensional models of massive gravity (see 
e.g. [SHE] for a discussion of this point). The inclusion of both the LCS term and the 
curvature-squared term of NMG leads to the "general massive gravity" (GMG) model [I]. 

Allowing for a cosmological term, and omitting an overall factor proportional to the 
inverse 3D Newton constant, the Lagrangian density for GMG, with spacetime metric 
9iiv (a* = 0, 1, 2), takes the form 



where the tensors G and S are, respectively, the Einstein and (3D) Schouten tensors, 
m and fi are two mass parameters, a is a dimensionless constant and A is the cosmo- 
logical parameter; in a maximally symmetric vacuum we have = —Ag^ u , where the 
cosmological constant A is given by 



Notice that Ao = allows a Minkowski vacuum, in which case perturbative unitarity 
requires a < 0; in other vacua it is useful to allow for arbitrary a (although for d/0 
one may rescale fields such that o 2 = 1). Linearization about the Minkowski vacuum 
leads to a generalization of the 3D Fierz-Pauli (FP) equations that propagates two spin-2 
modes of independent masses; NMG corresponds to the equal mass case, and TMG to the 
case in which one mode is decoupled by taking its mass to infinity. Linearization about 
other vacua leads to modified versions of these equations that depend on the ratios of 
the masses (/i, m) to the scale set by the cosmological constant. 

While a linear approximation usually allows a reliable count of the number of local 
degrees of freedom, there are cases in which it gives misleading results. A Hamiltonian 
formulation provides a way to count the number of local degrees of freedom without resort 
to linearization: this number may be defined as half the dimension per space point of 
the physical phase space (i.e. taking into account all constraints and gauge invariances) . 
For example, using the standard Arnowitt-Deser-Misner (ADM) Hamiltonian formalism 
for GR in a spacetime of dimension D, this definition tells us that there are D(D — 3)/2 
degrees of freedom, which coincides with the number of polarization states of a massless 
graviton. An extension of the ADM formalism to higher-derivative gravity theories in 
any spacetime dimension was worked out in [0] but it is not applicable to parity-odd 3D 




(1.1) 



A = -2m 2 a ± Va 2 + A 



A = A /m 2 . 



(1.2) 
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theories like TMG and GMG. The Hamiltonian formulation of TMG has been studied in 
many papers e.g. [TUHT2] but of most relevance here is the formulation that we shall refer 
to as the "minimal" formulation of Carlip [T3], some aspects of which were clarified by 
Blagojevic and Cvetkovic [14] via an implementation of Dirac's general procedure [16]. 

More recently, Dirac's procedure was used to find a Hamiltonian formulation of 
NMG [15], which was subsequently applied to the case in which A = —a 2 [17]. This 
case is special for many reasons [5] but the one of relevance here is that linearization 
yields equations with an accidental gauge invariance that describes a "partially mass- 
less" graviton with only one local degree of freedom [18j[l9]. The gauge invariance is a 
linearized Weyl invariance, which is "accidental" in the sense that it does not extend be- 
yond the linearized approximation, and for this reason one should suspect the reduction 
in the number of local degrees of freedom to be an artefact of the linear approximation. 
Indeed, it turns out that there is no reduction in the number of local degrees of freedom 
if this number is counted using the Hamiltonian formulation of NMG. This implies that 
NMG suffers from a linearization instability at "partially massless vacua". This result 
has been extended to GMG in [20J. 

There is one other case in which linearization leads to an accidental gauge invariance, 
which is again a linearized Weyl invariance. This is the case in which Ao = and 
a = 0; the linearized Weyl invariace is "accidental" because the non-linear theory is 
definitely not Weyl invariant [I]. The special subcase in which \x = oo was first analysed 
by Deser [2TJ, who showed that it propagates a single massless mode; we shall call 
this model "massless NMG". The generic case was analysed in [22], where it was called 
"topologically new massive gravity" (TNMG); this model was studied independently 
(under a different name) in [23]. Linearized TNMG propagates a single spin-2 mode of 
mass m 2 /fi, which becomes the massless mode of "massless NMG" in the \i — > oo limit. 
There is therefore an apparent reduction in the number of local degrees of freedom^ but, 
as in the "partially massless" case, we should expect this reduction to be an artefact 
of the linearized approximation. One of the aims of this paper is to use a Hamiltonian 
formulation to verify this. 

Implementation of Dirac's procedure for 3D massive gravity models leads to a rather 
large phase space with a correspondingly large number of constraints of both "first- 
class" and "second-class" in Dirac's terminology (first-class constraints being those that 
generate gauge transformations). Another purpose of this paper is to present a simple 
Hamiltonian formalism of NMG and GMG that, like Carlip's Hamiltonian formulation 
of TMG [13] . is "minimal" in the sense that (i) the only first-class constraints are the 
six that generate 3D diffeomorphisms and local Lorentz transformations, and (ii) the 
number of second-class constraints is minimized. 

Our starting point is a new action for GMG, and its TMG limit, that is "Chern- 
Simons-like" in the sense that it is the integral of a Lagrangian 3-form constructed by 
taking exterior products of independent one-form fields and their exterior derivatives, and 
so does not require a metric for its construction. As in Chern-Simons (CS) gravity models, 

1 The natural extension of TNMG to adS vacua is "chiral GMG", which occurs for A = —3a 2 + 
2p(p — 3a) — 2 (p — 2a) \/ {p — 2a) p, where p — m 2 //i 2 ; this is the GMG analog of the perturbative 
unitary limit A = — 3cr 2 of NMG [S]. Although there are discontinuities in the spectrum of linearized 
GMG at the chiral point, the linearized approximation does not lead to any reduction in the number of 
local degrees of freedom, except in the Minkowski limit where it degenerates to TNMG. 
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the one- form fields include the Lorentz frame, or dreibein, one-forms e a (a = 0, 1, 2) which 
can be used to construct a metric if we assume that its matrix of coefficient functions 
e M a is invertible. In a decoupling limit in which all propagating modes become infinitely 
massive, the action reduces to the Einstein-Cartan (EC) formulation of 3D GR, which is 
a CS gravity model p4,25j. In the partial decoupling limit that yields TMG, the action 
reduces to the sum of the EC action and the CS action for 3D conformal gravity of Home 
and Witten [25]. 

It might seem superfluous to say that the Hamiltonian formulation of massive gravity 
models breaks the manifest spacetime diffeomorphism invariance of the covariant action 
because the Hamiltonian formulation involves a distinction between time and space. 
However, the Hamiltonian form of the action for a CS theory is just the covariant action 
rewritten in a non-covariant way by performing a time/space decomposition of all fields. 
In contrast, the Hamiltonian form of the action of a CS-like massive gravity model cannot 
be obtained in this way; there are necessarily additional, secondary-constraint, terms 
that break this manifest spacetime covariance. While one such constraint is required for 
TMG [13], two are required for GMG. 

Using our Hamiltonian formulation of GMG, we are able to compute the dimension 
of the physical phase space and to see how this depends on the parameters of the model. 
In particular, we are able to show that the number of local degrees of freedom changes 
only in a decoupling limit and is therefore not discontinuous in any of the special limits 
discussed above, except in the limit that yields conformal 3D gravity. This means, in 
particular, that both the "massless NMG" limit of NMG and the TNMG limit of GMG 
suffer from linearization instabilities. 

2 Preliminaries 

As our covariant starting point for a Hamiltonian formulation of massive 3D gravity 
models is a Chern-Simons-like action for these models, we review here the aspects of 
the Chern-Simons gravity models that will be of relevance to us. This will also serve to 
introduce our conventions and terminology. 

2.1 Einstein-Cartan 

In the Einstein-Cartan (EC) formulation of 3D gravity, the independent fields are the 
dreibein e a (a = 0,1,2), a Lorentz-vector valued one-form, and an adjoint- valued local 
Lorentz connection one-form which we may trade for the vector- valued one- form u a . 
From these one may construct the Lorentz-covariant torsion and curvature two-forms 

T a = De a = de a + e abc u b e c , R a = du a + ^e abc u b u c , (2.1) 

where e abc is the antisymmetric invariant tensor of SO(l,2), with e 012 = 1. Here, and 
henceforth, products of forms should be understood as exterior products. The torsion 
and curvature two-forms satisfy the Bianchi identities 

DT a ee e abc R b e c , DR a = . (2.2) 
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This formalism allows the 3D Einstein-Hilbert action, with cosmological term, to be 
written in a first-order form as an integral of the Einstein-Cartan (EC) Lagrangian 3- 
form 

L EC = -ae a R a + ^e abc e a e b e c . (2.3) 
o 

As explained in the introduction, the constant a is included for later convenience; a = 1 
is the standard choice in 3D GR for our "mostly plus" metric signature convention. The 
EC action is already first-order in time derivatives and so constitutes a simple starting 
point for the Hamiltonian formulation of GR |27j . 

The field equation for uj a is T a = 0. If we assume that the dreibein matrix e p a is 
invertible, with inverse e a M , then we can solve for u a : 



-e-'e^ U,<V - -e/e^ d p e b = w/(e) , (e = det e/), (2.4) 



where e^ up is the totally antisymmetric invariant tensor density. Let us note here for 
future use that 

Under the same assumption of invertible dreibein, we have that 

l -e^ p R vp a = -ee u a G^, (2.6) 

where G pv is the Einstein tensor constructed from the metric := e p • e u = e p a e v h r\ab 
and the affine connection 

r M , A = co pa b e u a e b x + e a %e u a , (2.7) 

which becomes the standard Levi-Civita connection for the metric g pu when u = u(e). 
Using the fact that 2g ll ' J G flu = —R in 3D, where R is the Ricci scalar, we see that 
elimination of u a yields the Lagrangian density 

J2? = l -e [cR - 2A ] . (2.8) 

The identity (12. 6 j) may also be used to show that the e a -equation is, given u a = u) a (e), 
the Einstein equation aG^ = —Aog^, from which we see that Ao/cr is the cosmological 
constant for the EC model . 

A variant of the EC action, which will be useful for later purposes, is defined by the 
Lagrangian 3-form 

L EC = L EC + h a T\ (2.9) 

where the new Lorentz-vector-valued one-form h a is a Lagrange multiplier for the con- 
straint T a = 0. The equivalence of this action to the standard one can be seen by 
rewriting it in terms of e a and Q a = u a — h a jo\ one finds that 

Lec = L EC {e, n) + ^-e abc e a h b h c . (2.10) 

AO 

In this action the field h a can be trivially eliminated (given invertibility of the dreibein) 
whereupon the action reduces to the standard EC action. 
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2.2 Horne-Witten 



Now consider the Lagrangian 3-form 

Lhw = \oo a du a + \e abc u a u b uj c + h a (T a - be a ) + \bdb . (2.11) 
l o 2 

This is the CS theory for the 3D conformal group constructed by Home and Witten 
|26j . In addition to the one-form gauge potentials (e a ,u a ) of the EC model we have an 
additional Lorentz-vector valued one-form potential h a associated to proper conformal 
gauge transformations, and a Lorentz-scalar one-form potential b associated to local scale 
(Weyl) transformations. The relative coefficients are fixed, up to field redefinitions, by 
the requirement of local SO(2, 2) gauge invariance. 

The (h a ,u a ) fields are auxiliary fields in the HW model in the sense that they may 
be eliminated by their field equations, which are jointly equivalent to 

T a = be a , R a + e abc e b h c = . (2. 12) 

Assuming invertibility of e^ a , these equations imply that 

< = <(e, b) = <(e) + e™b v , V = "<We, 6) , (2-13) 

where S^ u (e, b) is the Schouten tensor for the connection o; At °(e, b); recall that, in 3D, 

Sfu, = Rpv — -g^yR, (2-14) 

where R^ is the Ricci tensor and R the Ricci scalar. Substituting for u using the first 
of equations f !2. 13j) . one finds that 

S^(e,b) = S^ + D^+^K-^b 2 ) , (2.15) 

where S^e) is the (symmetric) Schouten tensor for the standard torsion- free connection, 
and D indicates a covariant derivative with respect to this connection. Observe that 

S [ftv] (e,b) = d^b v ] , (2.16) 

from which we see that the second of equations (12.131) implies that 

fr[H = -d\p,b u ] . (2.17) 

This is also the 6-equation of motion. Finally, the e a -equation is 

(D + b)h a = 0. (2.18) 

Upon substituting the expressions for (h a ,co a ), one finds that all 6-dependence cancels 
and we are left with the equation C^ v = 0, where the left hand side is the symmetric and 
traceless Cotton tensor 

C^ = e^D T S p „(e). (2.19) 

This is the standard field equation of 3D conformal gravity. Off-shell equivalence of the 
HW model to conformal gravity can also be shown by using fl 2 . 1 3 [) to eliminate (h a ,u a ) 
from the HW action. One finds that the b field drops out and the resulting action is the 
Lorentz-Chern-Simons term for the composite connection u(e). 
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3 Massive 3D gravity 



Consider the following CS-like Lagrangian 3-form 

L G mg = -oe a R a + \k e abc e a e b e c + h a T a + ' 

o 



2/i 



u a du a + -e abc u a u b u c 



f a R a + ^ Cafbfc 



(3.1) 



In the limit in which both /i — > oo and m — > oo we get the variant EC Lagrangian 3-form 
of (12.91) . In the limit that m — > oo for finite /i we get the Lagrangian 3-form 



j tmg = -ve a R a + \k e abc e a e b e c + h a T a + ^~ 
6 2/x 



u a dco a + ^e abc co a UJ b UJ c 



(3.2) 



which is known to describe TMG [T21IT3] . 

In the generic case of finite fi and finite m the above Lagrangian 3-form is equivalent 
to the one proposed for GMG in [T5|[2"0] if one assumes an invertible dreibein field, but it 
seems not to have been previously appreciated that the GMG Lagrangian can be written 
in the above way. Let us first confirm that it does indeed describe GMG. The (h a , u a , f a ) 
equations are jointly equivalent to 



T a = 



which imply that 



and hence that 



~Df a + -R a + e abc e b h c = , R a + e abc e b f c = 



77T 



h 



[H 







(3.3) 

(3.4) 
(3.5) 



Back substitution yields the Lagrangian density 



1 



GMG 



1 



e laR - 2A + G^{e)S„ v {e) + -JSfj 



1 



'LCS 



(3.6) 



where J^lcs is the standard Lorentz-Chern-Simons term of 3D conformal gravity. For 
a = — 1 this defines GMG, and NMG is obtained in the /i — > oo limit. 



3.1 Hamiltonian preliminaries 

The action defined by integration of (13. ip over a 3-manifold with a Cauchy 2-surface is a 
starting point for a Hamiltonian formulation of GMG, and its various limits. Observe that 
the Lagrangian density corresponding to the Lagrangian 3-form (13. ip takes the general 
form 

1 i 

& = ^9r S a r ■ da 3 + -f rst a r ■ a s x a 1 , (3.7) 
z o 
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where a r (r = 1, . . . , N) are N 'flavours' of 3- vector one-forms, and we use a 3D Lorentz- 
vector algebra notation in which, e.g. 



a P ■ a q = V ab <a q b , (a s x a 1 ) a = e abc a s b al , (3.8) 

where r\ is the Minkowski metric (of "mostly-plus" signature) and e is the antisymmetric 
Lorentz-invariant tensor such that e 012 = 1. The constant coefficients g rs , f rs t are defined 
such that they are symmetric under exchange of any two indices. We can view g rs as a 
metric on the 'flavour' space if we assume (as is the case for the models of interest) that 
it is invertible. We need the N = 3 case for TMG and the N = 4 case for GMG. The 
time/space decomposition 

a r = dta r + dCa r i (* = 1,2) (3-9) 
yields a Lagrangian density of the form 

SC = -~r) ab evg rs a\ a d° h + a r J a r , (3.10) 

where 5 y = e 0iJ . The phase space has dimension QN. The time components a r impose 
3N (primary) constraints <p r = 0. 

At this point we can anticipate how the Hamiltonian formulation will lead to the 
conclusion that GMG describes two local degrees of freedom, realized in the linear theory 
as two massive spin-2 modes. We can see from ( 13. 5ft that the equations = and 
f\ij] = are additional secondary constraints, since they are conditions on canonical 
variables that are not imposed by Lagrange multipliers but which are a consequence of 
the equations of motion. Following Carlip's analysis of TMG [13] . we are led to consider 
the modified Lagrangian 

&GMG+ = &GMG + bo e^hij + c f tJ (3.11) 

The phase space is spanned by the 2-space components of the four 3- vector one-form fields 
and hence has dimension per space point of 4 x 3 x 2 = 24. The time components of the 
fields impose 4 x 3 = 12 primary constraints, to which we must add the two secondary 
constraints imposed by the new variables (&o,Co), making a total of 14 constraints. Of 
these we expect six to be first class, corresponding to the built-in local Lorentz and 3- 
space diffeomorphism invariance. The physical phase space will then have dimension, 
per space point, of 24 — 14 — 6 = 4, corresponding (by our earlier definition) to two local 
degrees of freedom. 

We shall fill in the details later. First we wish to address two other important ques- 
tions: 

1. Is the modified Lagrangian, with the extra secondary constraints, equivalent to the 
original unmodified Lagrangian? 

2. Is the modified Lagrangian still invariant under 3-space diffeomorphisms? Although 
these terms respect manifest 2-space diffeomorphism invariance they break manifest 
3-space diffeomorphism invariance. 
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Notice that we could restore manifest 3-space diffeomorphism invariance by starting 
from the Lagrangian 3-form Lqmg — be-h — ce-f, so that the secondary constraints are 
imposed by the time components of the new one-form fields b and c. There is still the 
issue of equivalence of the new field equations to the original GMG field equations; in 
fact, equivalence is lost in making this modification but it can be restored (subject to a 
reservation to be discussed below) by adding kinetic terms for (b, c) to arrive at the new 
Lagrangian 3-form 



J GMG++ 



-ae a R a + \k^ abc e a e b e c + h a (T a - be a ) + -Ucft +±- c {db-e- f) 
6 2ii m z 



+ 



2li 



1 

— t 

3 



bj a duj a + ^-e abc u a u b u c 



77T 



faR a + -e abc e a f b f c 



(3.12) 



(3.13) 

, c) fields yields the minimal 



Observe that 

^GMG++\b l =c i =0 = ^GMG+ ■ 

In other words, setting to zero the space components of the 
modification of ( 13. lip . 

Whether we choose to consider the GMG + or the GMG ++ modification of the GMG 
action, we still need to address the issue of whether the modified field equations are 
equivalent to those of GMG. This is not obvious even for TMG; it was observed in [13] 
that Carlip's Hamiltonian formulation of TMG (which amounts to using TMG + ) was 
incomplete because it did not include a proof of equivalence to TMG. We provide a 
proof here that the field equations of both TMG + and TMG ++ are equivalent to those 
of TMG. As we shall see, the analogous equivalence proof for GMG is more involved 
and does not apply without qualification. We shall begin our analysis by considering 
the 3-space covariant "++" modifications since this is simpler and the analysis is easily 
adapted to the non-covariant "+" modifications. 



TMG ++ . By taking m — > oo in (I3.12p we arrive at the following Lagrangian 3-form 
in terms of the one-form fields of the HW model 



L 



-ae a R a + lA e abc e a e b e c + h a (T a - be a ) + -Ud& 
6 2il 



1 

+ 2^ 



u a duj a + -e abc u a uj b uj c 



(3.14) 



In the limit that o — > we must also set Ao = in order to get consistent field 
equations, and in this case we recover the HW model after a rescaling of h a . In 
other words, the model under consideration is essentially defined by an action that 
is the sum of the EC and HW actions. We shall now show that this is another 
description of TMG. 

Consider first the (h, u) equations, which are jointly equivalent to 



T a = be a , 
These may be solved to give 

< = <(e,6) 



-R a + e abc e b h c - abe a 







h 



,„> - Spvfa b) + (je^ vX b X ■ 



(3.15) 



(3.16) 



S 



The second of these equations implies that 

fr[H + -d^K] = (re^ uX b x . (3.17) 
/i 

However, the 6 equation is equivalent to the vanishing of the left hand side, from 
which we deduce (assuming a ^ 0) that 6 = 0. Therefore, the combined (h a ,u a ,b) 
field equations can be solved for these fields to give 

w/ = w/(e), V = -<V( e ), ^ = °- ( 3 - 18 ) 

Back-substitution then yields the standard TMG action. Alternatively, we may 
observe that since the field equations that follow from (I3.14p imply that 6 = 0, the 
remaining equations are equivalent to those of TMG; specifically, those that follow 
from the TMG limit of the GMG Lagrangian 3-form (13. ip . 

GMG ++ . We now turn to the general case of the Lagrangian 3-form (I3.12p . The 
(h, f) equations 

T a = be a , R a + e abc e b f c = ce a (3. 19) 

can be solved to give 

= <(e, b) , = -S^e, b) + e^cf , (3.20) 

and hence 

/[H + d [fl b v] = -e^ p c p . (3.21) 

On the other hand, the c-equation implies that f\p U ] + d[^b u ] = 0, and hence that 
c = 0. 

Let us now consider the w-equation; this is equivalent, given (I3.19P and c = 0, to 

_L Df a _ e -bc eb U _ I A + abe a = g ( 3 22) 

m l \ fi J 

which has the following solution^: 

m 2 ^V-^/^) = ^( e ) + ^(e P<Tr ^A) + e/ A 6 T ^(e) 

+ t v rX b T D x bp + i e ^ A 6 A (b 2 - 2am 2 ) . (3.23) 

However, the b and c equations combined tells us, given c = 0, that the antisym- 
metric part of the left hand side is zero. So the antisymmetric part of the right 
hand side is zero, and this is equivalent to the equation 

~} w b v = , E^ = GT - m v - (2am 2 + D ■ b - b 2 ) . (3.24) 

Although 6 = solves this equation, there are other possible solutions, so GMG ++ 
is not strictly equivalent to GMG. Instead, it appears that GMG is equivalent to 
one "branch" of the model defined by the GMG ++ action. 



2 It is useful to use here the fact that (D + b)f a is 6-independent given (|3.19[) . To see this, compare 
with the equation of motion (|2.18|) of the HW model. 
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This branch equivalence goes beyond the statement that the field equations of 
GMG ++ reduce to those of GMG when we choose the 6 = solution of (I3.24p . An 
analysis of fluctuations about any solution of GMG ++ with 6 = leads to the linear 
equation [G^ IV — 2am 2 g fiu ] 5b u = 0, which implies that 6b = 0, generically, and hence 
that a linear stabilty analysis for solutions of GMG ++ with b = is equivalent to 
a linear stability analysis in the context of GMG. An exception to this state of 
affairs occurs when we consider fluctuations about a maximally-symmetric vacuum 
with G^ v = 2am 2 g^ u . This is precisely the case in which the massive gravitons 
become "partially massless" [I]. In this case b is undetermined at the linear level, 
so there must be an "accidental" gauge invariance of the linear theory that allows 
it to be "gauged away" . At the non-linear level it is still present but we may also 
still choose the 6 = solution. 



Recall that we introduced the GMG ++ Lagrangian 3-form (13.121) in an attempt to 
restore the manifest 3-covariance that is lost when considering the minimal modification 
of GMG + . However, if we were to take this as our starting point for a Hamiltonian 
formulation we would need to include the non-covariant conditions 6j = q = as new 
secondary constraints because these are constraints on canonical variables that are either 
implied by the field equations or imposed consistently with them. As a check we observe 
that we would then have a phase space of dimension 4x3x2 + 2x2 = 28, per space 
point, and 12 + 2 + 4 = 18 constraints. Given that six of these constraints are first class 
we then get a physical phase space of dimension 28 — 18 — 6 = 4, as expected. However, 
this would give us a non-minimal Hamiltonian formulation that has no obvious advantage 
over the minimal formulation provided by GMG + . 

We therefore return to the GMG + Lagrangian density of (13. lip . Although the in- 
clusion of the secondary constraint terms, with Lagrange multipliers 6o and Cq, breaks 
manifest 3-space covariance, the field equations will still be 3-space covariant if they can 
be shown to be equivalent to those of GMG, so we now need to address this issue, which 
we can do by adapting our analysis above for GMG ++ . We shall again consider the TMG 
case separately. 

• TMG + . As bi = the 6 equation is now huj] = 0. The (h a ,u a ) equations are just 
those of TMG ++ but with 6j = q = 0. This applies, in particular, to (I3.17P from 
which we see that ab is proportional to huj-i, and is therefore zero. As long as a 
is non-zero this implies that 6 = 0, and hence 6 = 0. Using this, the remaining 
equations become those of TMG. 

• GMG + . By setting 6j = q = in f )3.2ip we deduce that /[ij] is proportional (for 
finite m) to Co, but the cq equation is /[y] = 0, so we deduce that Co = 0, and hence 
that c = 0. So now we have equations that are equivalent to those of GMG except 
for possible 6o dependence. 

The combined (6o,co) equations imply the symmetry of fihij — /y. Using this in 
( I3.23P and setting 6j = we arrive at the equation 



where H is the tensor of (I3.24p . This equation has 6o = as a solution, and choosing 
this solution we get equivalence with GMG. 
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We see that the branching of possibilities for b that we found in our previous analysis 
of GMG ++ has not been entirely eliminated. The equation (I3.25P has 60 = as one 
solution but there are other possibilities, at least one of which coincides with the b Q = 
solution at "partially massless" vacua. A similar branching of possibilities at these vacua 
was found in the analysis of [20]. From our perspective, it appears necessary to insist on 
60 = because otherwise the field equations are not those of GMG. 



4 Hamiltonian formulation 



Following the discussion of the previous section, we take as our GMG Lagrangian 3-form 



L G mg+ = -ae a R a + \k e ahc e a e b e c + h a T a + ^~ 

o 2/i 



u a du a + -e abc u a u b u c 



1 



77T 



f a R a + -e abc e a f b f c 



be ■ h -ce ■ f . 



where 



b = b dt , c = codt . 
The (bo, Co) fields are now Lagrange multipliers for the constraints 

h m = ' fm = • 



(4.1) 
(4.2) 
(4.3) 



We shall call these the "secondary" constraints since they are secondary in the context 
of the GMG action fl3~T|) . 

The Lagrangian density still takes the general form (13. 7p except that we must now 
add the secondary constraints. After a time/space decomposition we then arrive the 
Lagrangian density 

Sf+ = --e^grs^ ■ a) + a r ■ <p r + 6^1 > ( 4 - 4 ) 
with I — 1, . . . n, and n — 1 for TMG and n = 2 for GMG. The constraint functions are 



e 13 dia)g rs + -/ rjrt oJ x a] 



3 / ' 



^ = —fj iPq A pq 



where fi tPq = —fi, qp is a new set of constant coefficients, and 



A pq = e ij a p ■ a q - . 

1 J 



(4.5) 



(4.6) 



Observe that A pq = —A qp . The quadratic term of (14.41) gives us the Poisson brackets 

{< a (t;),a* b (0} pB = r ]a ,g rs e ij 6W(t;-() , (4.7) 

which we may use to compute the matrix of Poisson brackets of the constraint functions. 
To this end, it is convenient to first define 



(4- 
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where the test functions a r a are arbitrary except that we choose them such that no surface 
terms arise upon integration by parts. A calculation using (14. 7p yields the result 



{0(a),0(/3)} PB = 0([a,/3]) + A[r/.W d 2 £a r -P a A™ 

+ 2f r[s f q]pt [tPtaffiiV*)" , (4.9) 



where 

[a, 0} c t = e abc a r Mst , Vg = e^a\ a a) h . (4.10) 

As we discuss in the following subsections, it turns out that for the models of interest 
the A terms in these Poisson brackets are all zero as a consequence of the secondary 
constraints. It follows that, on the constraint surface, 

W(e),rf(0} = ^-o, (4.11) 

where 

P?.{t - = frisUpt (V ab ) m 5 {2) (£ - C) • (4.12) 

This 3N x 3N matrix P plays a crucial role in the analysis to follow. However, we will 
also need to take into account the Poisson brackets of the primary with the secondary 
constraints. A calculation for the generic model shows that 

{0(a), V/}ps = e ij [dioralfm - «X * a'/r///^] • (4-13) 

For NMG and GMG we also need the Poisson bracket of the two secondary constraints; 
we shall see that these two constraints are in involution. 



4.1 TMG 

Separating the quadratic from the cubic terms in (13. 2p . we find that 
Ltmg = -<ye a du a + h a de a + ^-u a du a 

+ e abc \h a oo b e c - -e a co b Lo c + —oo a uj b uj c + — e a e b e c \ . (4.14) 
L 2 6[i 6 J 

We can simplify the quadratic term by setting 

co a = Q a + afie a , h a = k a + V/ie". (4.15) 
In terms of (e, k, Q) the Lagrangian 3-form is 



Ltmg — k a de a + —Q a dQ a + e ahc 
2/x 



k a Q b e c + —Q a Q b Q c + an k a e b e c + -A e a e b e c , 
6/x 6 

(4.16) 

where 

A = A + a 3 /? . (4.17) 
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Apart from differences in notation, this result differs from the analogous result of the 
HW model only in the a- and A -dependent cubic terms. 

We are now in a position to write down a Hamiltonian form of the action, by per- 
forming a time/space split in ( I4.16P and adding the secondary constraint. This gives 
us, 

^tmg+ = • ej + ■ Uj | + e • e + ujq ■ W + k ■ (p k + b A ek , (4.18) 



where 



(j) e = e 1 -' <^Dikj + 2afiki x ej + -A ti x ej 

4>k = ^ + cr^ei x ej} . (4.19) 

These constraint functions are, of course, just the specialization to the case in hand of 
those given by the general formula (14. 5p . We see from this result that the phase space, 
spanned by the space components of the Lorentz 3- vectors (cu, e, k) has dimension 18 per 
space point, and that there are 10 constraints. 

We actually have no need for the above explicit expressions for the primary constraint 
functions r since we may use the general result for their Poisson brackets that we have 
already computed in terms of the various coefficients that define the model. To do this 
we first read off from (I4.16P the non-zero components of g and /, which are 

9ek = 1 , 9nn = At" 1 , (4.20) 

and 

fkde = 1 , fnnn = ^ 1 , fkee = ^crfx, f eee = A , (4-21) 
from which it follows that 

f Q ke = ^ , f Q nn = f e ne = f k nk = 1 , 

fee = A , f\ e = f ke = 2afi. (4.22) 

Using these results we find that the 9x9 P-matrix of ( 14.111) takes the following form 
in the (Q, k, e) basis: 

<«.)„ = ( o Q ) <«■») 



where Q is the antisymmetric 6x6 matrix [f| 



V'ee T/efc 



Q = /^ (2) (£-C) _«\ k • (4.24) 

V V ab V ab J 

Note that this matrix is independent of both a and Ao- The zeros of the first row 
and column of P (actually 3 rows and 3 columns because we suppress Lorentz indices) 



3 It is antisymmetric because — (V^) T = V^jf. 
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are expected from the built-in local Lorentz invariance, which ensures that the Poisson 
bracket of 0q with any other constraint is zero on the constraint surface, i.e. that the 
constraints 0^ are first-class. We may use this built-in local Lorentz invariance to choose 
a local frame for which e\ = (010), = (0 01), in which case the secondary constraint 
implies that k\ = k\. It can then be easily verified using Mathematica that the matrix 
Q (and hence P) has rank 2. 

However, we still have to take into account the one secondary constraint, with con- 
straint function if) = A ek . From the general result ( 14.13P we find, in this instance, that 

{0(a), 4>} PB = e ij (Diafkj - Diofej) + (2afia h + A a e ) e^'e* x e 5 (4.25) 

where we use the shorthand 

Diofctj = dia r a S j - a r Vti x a* . (4-26) 

The absence of any term involving a Q is expected from the fact that the secondary 
constraint function is a Lorentz scalar, but both e and 0^ have non-zero Poisson brackets 
with if): 

{MZ)MC)}p B = -e y {W (2) (£-0- [fiiX^-Ao^xejJ^e-C) 

{MZ)MC)}p B = e ij {e j d i 5^(C-C) + [^iXe j + 2a f ie i xe j }5 ( - 2 \C-0} ■ (4.27) 
This shows that the 10 x 10 matrix P of Poisson brackets of constraints takes the form 

P = ( n IV (4-28) 



n 

where Q is a 7 x 7 antisymmetric matrix of the form 

Q v \ _ ( {4>k,ip} 



V J ~ \ {0e,^} 



l PB 

^PB 



(4.29) 



All dependence on a and A enters through the column vector v. If this column vector 
is in the column space of Q then the rank of Q equals the rank of Q, i.e 2. This would 
imply that there are 8 first-class constraints and hence 18 — 10 — 8 = local degrees of 
freedom. This special case is realized if and only if a = Aq = 0, as expected because this 
is the limit in which TMG degenerates to conformal gravity, which has no local degrees 
of freedom. The 8 gauge invariances are what remains of the local conformal invariance 
in the gauge in which 6, = 0. In all other cases v is not in the column space of Q, and 
the rank of Q is then 2 greater than the rank of Q, i.e. 4. This implies that there are six 
first-class constraints, per space point, corresponding to six gauge invariances that can 
be identified as those of 3-space diffeomorphisms and local Lorentz invariance [13]. The 
dimension, per space point, of the physical phase space is therefore 18 — 10 — 6 = 2, as 
expected for TMG. 

4.2 NMG 

For a reason that will become clear, we develop the Hamiltonian formalism for NMG 
separately from that of GMG. Taking the \i — Y oo limit in (13. ip we arrive at the NMG 
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Lagrangian 3-form: 
Lnmg — 



-ae a R a + lA e abc e a e b e c + h a T a - i- 



6 



f a R a + -e abc e a f b f c 



. (4.30) 



Separating the quadratic and cubic terms, and then simplifying the former by defining 
the new variable 

-3/*-^, (4-31) 



7T 



we arrive at the Lagrangian 3-form 

Lnmg = Kde a + n a duj a + e abc 



h a u b e c + -ir a u°u c 
2 



a, ,b, ,c 



+ ^abc 



Tfl l y, 

e a n b 7r c - m 2 ae a e b % c + -A - e a eV 

2 6 



(4.32) 



where 



A = A - 3m 2 a . (4.33) 

Making a time/space split, and then adding the two secondary constraints, we arrive at 
the NMG Lagrangian density in Hamiltonian form 

&nmg+ = -e^ihi-e.+^-io.j + boA^ + coA^ 

+ uq ■ 4>u + e • <p e + h ■ 4> h + ttq ■ 07T , (4.34) 



where 



6 lj {DiTCj + d X hj} 



3 i ' 



>) h = e l W iej , 



Tfl 1 " 

e tJ ^ Dihj — x n j — 2m 2 a x ttj + —A e, x 



DiUJj + -a;, x 



— m 2 ei x 7Tj — mVej x 



(4.35) 



The phase space now has dimension, per space point of 24 but there are a total of 14 
constraints. Our next task is to determine the dimension of the subspace of first-class 
constraints and hence the number of gauge invariances. 

As for TMG, we do not need to use directly the above explicit expressions for the 
primary constraint functions because we may instead use the general result for the Poisson 
brackets of the constraint functions that we computed earlier. For this we need the 
expressions for the coefficients g and / in the (u, it, h, e) basis, which we may read off 
from (14.321) . The non-zero components of g and / are 



eh 



9t, 



fhuie fir 



u 



from which it follows that 



-m 



f, 



-2m 2 a . 



fe 



A 



(4.36) 



re rh ruJ ptt en i rh £U 

J uje J uh J ujuj J eh J ttlj j J tttv J 

fen = r ee = -2mV, f k ee = A . 



-m 



(4.37) 
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Using these results we find that the P matrix of (14.111) in the (u, tt, h, e) basis again takes 
the form (14.231) but now with a 9 x 9 antisymmetric submatrix Q; suppressing Lorentz 
indices, we have 






yee 


_yeh 


yee 





_yen 


_yhe 


_yne 


yhn _|_ y 



Q = mH^\i-Q I V ee -V en ]. (4.38) 

Note that this matrix is independent of both a and Ao. It is antisymmetric because, for 
example, the transpose of V ee is — V ee . For the same choice of frame for ei a that we used 
for TMG, we have h\ = h\ and it\ = tt^, and a Mathematica calculation can then be 
used to show that Q has rank 4. 

Now we must take into account the two secondary constraints, with constraint func- 
tions 

ih = A eh , ^2 = A ew (4.39) 

It is easily verified that {^i, 1^2} pb ^ ^2, so this PB is zero on the constraint surface. 
The non-trivial Poisson brackets are 

{(p(a),ipi}p B = e 1 ^ [Dia e hj — D^ej — vr^oFiii x ej] 



+ £ 



'■J 



(a e + a n ) A ei x ej — 2m a a e Hi x ej 

{0(a), ^2} PB = e ij [Aa e 7Tj - D^ej + a e hi x ej + a h d x ej (4.40) 
where we again use a shorthand notation: 

Did' a) = dioTa* - a r Ui x a* . (4.41) 



This gives us a 14 x 14 P- matrix of the general form (I4.28P but now with an 11 x 11 
antisymmetric sub-matrix Q of the form 

/ Q Vi V 2 \ / {(pn,^l}pB \ 

Q= -vf , vj=\ {Mj} pb , (4.42) 

V ~ V l / V iteMpB ) 

where the Poisson brackets {4>ri'4 , i}pB can read off from (14 .40 p . Observe that the 
components of these column vectors are sums of terms that are either linear or quadratic 
in canonical variables, and that all dependence on a and Ao is contained in the quadratic 
terms. 

For a matrix Q of the above form, its rank equals the rank of Q (i.e. 4) if both v\ 
and v 2 are in the column space of Q, and it equals the rank of Q plus 4 (i.e. 8) if V\ and 
v 2 are linearly independent and no linear combination of them is in the column space of 
Q. In all other cases the rank of Q is the rank of Q plus 2 (i.e. 6). If the quadratic terms 
of vj were absent, then both these vectors would be in the column space of Q, so that 
Q would have rank 4. However, the most we can do to eliminate these quadratic terms 
is to set ex = Ao = and this still leaves some quadratic terms, which are sufficient to 
ensure that the column vectors Vj are both independent and that no linear combination 
of them is in the column space of Q, so the rank of Q is 8 independently of the values of 
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a or A . This means that 6 of the 14 constraints are first-class, as expected, and hence 
that the dimension of the physical phase space, per space point, is 24 — 14 — 6 = 4. This 
is consistent with the linearized analysis of the generic NMG model, which shows that 
there are two propagating modes, and with the Hamiltonian results of [15], but it also 
applies in the a = Aq = limit that yields "massless NMG", and in that case it is not 
consistent with the linearized analysis of [21]. We conclude that "massless NMG" suffers 
from a linearization instability. 

4.3 GMG 

For GMG we proceed initially as for NMG, separating the quadratic from the cubic terms 
of the GMG Lagrangian 3-form (13. II) and then making the change of variable ( I4.3ip to 
get to 



^gmg = h a de a + ^-u a du a + 7r a du a + e abc 

2yU 



La, ,b_c i . ,a, ,b, ,c i ^ _a. ,b, ,c 

n u e H u cu cu H — n u uj 

6fi 2 



77? 1 ~ 

e a TT b 7T c - m 2 ae a e V + -A e a e b e c 

2 6 



(4.43) 



To further simplify the quadratic term we now set 

, ,a r\a A* a 
OJ =16 -71 , 

m 2 

where fl a is a new independent connection. This gives us the Lagrangian 3-form 
Lgmg 



(4.44) 



h a de a + —Q a dQ a - --n a dix a + e abc 
2/i 2 



h a Q b e a + —Q a Q b Q c - ^7T a Q b TT c 
6/i 2 



+ tabc 



U Tfl 1 ^ 

— vr a 7rV - /i/ i a 7r fe e c e a 7r b 7r c - m 2 ae a e b n c + -A e a e b e c 

3 P 2 6 



(4.45) 



Observe that all f2-terms in the cubic term covariantize the quadratic terms with respect 
to local Lorentz transformations, so local Lorentz invariance is still manifest. The 'flavour' 
space metric is simple in the new (Q, tt, e, h) basis, but it is no longer simple to consider 
the /i — > oo limit that yields NMG. This is why we first dealt separately with the NMG 
case; having done so we may now assume that /i is finite. 

Making a time/space split and adding the two secondary constraints, we arrive at the 
following Hamiltonian form of the GMG Lagrangian density: 



GMG+ 



s tJ \ hi ■ e + — a ■ Q 



b A eh - c a A e 



+ Uq ■ W + e • e + h ■ <fi h + 7T 



(4.46) 
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diflj + x Qj 



7 ^ 

+ Ci X hj - -7Tj X 7Tj . 



e u < Dj/ij — fj,hi x TTj 7Tj x 7tj — 2m 2 a x 7Tj H — A x 



where 

h = e tj {Diej - fi-Ki x ej} , 

m 2 2 I 

2 J 2 

0^ = | — ^DiTCj + /iV.; x 7Tj — m 2 e.i x 7Tj — m 2 aei x e-,- j . (4.47) 

Again, we do not need to use these expressions directly because we may instead use 
the general result for the Poisson brackets of the constraint functions that we computed 
earlier in terms of the coefficients g and / that define the model. From (I4.45P we see that 
the non-zero coefficients in the (Q, e, tt, h) basis are 

Qeh = 1 , 9nn = /i" 1 , 9nn = -fJ> , (4.48) 
and that the non-zero components of f rs t are 

fhQe — 1 ) fnnn = ^ 1 , fnnfl = — A* , 

fnnn 2/i , fhire A* j 

fenn = -m 2 , / ee7r = -2mV , / eee = A . (4.49) 

It follows that 



f e Cle 


— J nh 


P7T 

— J flTT ~ 


J 


= 1, 


J he — A 1 , 


f n - -u 2 

Jnn h 1 i 


rn 

J he 


= 1, 


T7T 

J 7T7T 


-2/i, 


fn 

J en 




f\ e = 2mV//i , 


re 

J ne 


- f h 

J nh 


= -//, 


J 7T7T 


= —m 2 


J en 


-2m 2 a, 


fh 

J ee 


= A . 













(4.50) 

Using these results we find that the P matrix of (14.111) in the (Q, tt, h, e) basis again takes 
the form (14.231) but the 9x9 antisymmetric submatrix Q is now 



2 

m 



m 2 V ee fxV ee -fiV eh - m 2 V ew 

Q = — <5 (2) (£-C) I ^V ee -nV e7T ]. (4.51) 

V \ -^y he - m 2 V we -fxV ne fx(V wh + V h7T ) + m 2 V nw 

Once again, this matrix is independent of both a and Ao- Using the same choices for 
as before, we can now use Mathematica to evaluate the rank of this matrix. We assume 
that neither m 2 nor fi is zero or infinity. The result is that Q has rank 4. 

Now we consider the secondary constraints; the constraint functions are again 

^ 1 = A e \ ij 2 = A e \ (4.52) 
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and it is again straightforward to verify that {V'ijV^Ipb oc i/) 2 , so this PB is is zero on 
the constraint surface. Next, we compute 



{0(a), *l>i} PB = 

{(p(a),i) 2 }p B = 
+ e ij 

where, as for TMG, 



r'J 



DiQ! e hj — Dia h ej + fia e iTi x hj + [Aa e — 2m 2 aa K ) x e,- 



+ e tj [fia n ei x hj — (m 2 ^ + fia h + 2m 2 aa e ) iTi x ej] , 
e lj [Dia e iTj — Dia n ej + fia e iTi x 7Tj + a e /ij x ej] 



(4.53) 



a e — /ia" J x 7Tj + ( a + 



2m 2 a 



m 



a H a' 

/i /i 



6 j X 



Diofaj = dia r a S j — a r Qi x a* 



(4.54) 



This again gives us a 14 x 14 P-matrix of the general form (14.281) with an 11 x 11 
antisymmetric sub- matrix Q that is of the general form ( 14.421) . but with column vectors 
vj that we now read off from (|4.53|) . 

Once again all dependence on both a and A is contained in those terms in vj that are 
quadratic in canonical variables. In the absence of these quadratic terms both vectors vj 
would be in the column space of Q and the rank of Q would then be the same as the rank 
of Q (i.e. 4). However there are quadratic terms, and however we choose a and Ao the 
vectors vi are linearly independent and no linear combination of them is in the column 
space of Q. The rank of Q is therefore 8, independently of the values of a or Aq. As 
for NMG, this implies that the dimension of the physical phase space, per space point, 
is 24 — 14 — 6 = 4. This is consistent with the linearized analysis of the generic GMG 
model, and with the Hamiltonian results of [15J . but it also applies in the a = A = 
limit that yields TNMG, and in that case it is not consistent with the linearized analysis 
of [221123] ■ We conclude that TNMG suffers from a linearization instability. 



5 Discussion 



We have shown that the action for the 3D "general massive gravity" (GMG) model 
incorporating both "topologically massive gravity" (TMG) and "new massive gravity" 
(NMG) [4] can be written as the integral of a Lagrangian 3-form constructed from one- 
form fields (including a 3-vector dreibein) and their exterior derivatives. The action is 
then defined without the need for a metric, or even a density. We should stress that 
this reformulation of 3D massive gravity models depends on the special combination of 
curvature-squared invariants that appear in the NMG/GMG action, which is remarkable 
because this combination was not invented for this purpose. It would be interesting to 
see if a similar metric-independent formalism is possible for 3D massive supergravity 

|22I|2SII22| . 

We have called these metric-independent actions "Chern-Simons-like" because of their 
similarity to Chern-Simons (CS) theories of gravity. Strictly speaking, they define a 
generalization of 3D massive gravity models because equivalence to the usual actions can 
be established only if the dreibein field is assumed to be invertible. This is also how 
Chern-Simons (CS) theories of gravity become equivalent to standard metric theories of 
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gravity. The difference, from the perspective of this paper, is that CS theories require 
special coefficients for the various terms in the action, with the result that there are no 
local degrees of freedom. 

The absence of local degrees of freedom in CS gravity models is also apparent from 
their Hamiltonian formulation, which can be found directly by a time/space decompo- 
sition; the phase space dimension, per space point, is exactly twice the number of local 
phase-space constraints, all of which are "first-class", so the dimension, per space point, 
of the physical phase space is zero. In contrast, additional constraints are needed for the 
Hamiltonian formulation of Chern-Simons-like models. A single additional, "secondary", 
constraint suffices for TMG [13j[Tl], and this leads to the what we have called the "min- 
imal" Hamiltonian form of this model. We have used the Chern-Simons-like formulation 
of NMG and GM G to find an analogously "minimal" Hamiltonian formulation requiring 
two additional constraints. 

Secondary constraints are needed in the Hamiltonian formulation whenever the field 
equations imply constraints on the canonical variables that are not already imposed by 
the time-components of the one-form fields used to construct the Lagrangian. Since these 
are constraints on the space-components only, imposing them via Lagrange multipliers 
leads to an action that is no longer manifestly invariant under 3-space diffeomorphisms 
in the sense that (in contrast to the CS case) it is not the time/space decomposition of a 
manifestly covariant Lagrangian. This is not a problem if the new field equations imply 
the vanishing of the new Lagrange multipliers, because the new field equations are then 
equivalent to the original ones. For TMG it is easy to show that this is precisely what 
happens. In the NMG/GMG case, however, the new Lagrange multipliers are zero only 
on one branch of the solution space of the new equations, with other branches meeting 
it at "partially massless" vacua, as has also been found in the Hamiltonian approach 
of [20] ■ Our investigations led us to conclude that equivalence to NMG/GMG holds only 
on this one branch. 

Fortunately, this "branch equivalence" of our Hamiltonian formulation of NMG and 
GMG is sufficient for our main purpose, which is a determination of the number of local 
degrees of freedom of these models. The results of this computation, for all of the models 
considered in this paper, are summarized in the following table, where each model is 
defined by the combination of invariants that it includes; these are the Einstein-Cartan 
(EC) action for 3D GR, the Lorentz-Chern-Simons (LCS) term of 3D conformal gravity 
and the NMG curvature squared-invariant which, by itself, defines "massless NMG" 
(mNMG): 



Lec 


Llcs 


LmNMG 


Name 


Degrees of Freedom 


X 






Einstein-Cartan 







X 




Conformal Gravity 





X 


X 




TMG 


1 


X 




X 


NMG 


2 






X 


massless NMG 


2 


X 


X 


X 


GMG 


2 




X 


X 


TNMG 


2 



20 



These results confirm those of p~5|!17|J20j . in particular the absence of any discontinuity 
in the number of local degrees of freedom in the special case in which a linearized analysis 
yields "partially massless" gravitons. There are other special cases, however, and we have 
focused on the limit of GMG that yields "topologically new massive gravity" (TNMG), 
where a linearized analysis also exhibits an apparent reduction in the number of local 
degrees of freedom. Our results show, as expected, that this is an artefact of the linearized 
approximation and hence that TNMG suffers from a linearization instability. 

Our result for TNMG also applies to its parity-preserving "massless NMG" limit. 
This model was argued in [21 J to be renormalizable but the argument depends on the 
accidental linearized invariance of the linearized theory. This gauge invariance is "acci- 
dental" because the non-linear theory is certainly not Weyl invariant, although it does 
have a "conformal covariance" property [4] that explains why the linearized theory is 
linearized Weyl-invariant [30]. Similar considerations explain why linearized TNMG also 
has an accidental linearized Weyl invariance. 

From this discussion, it should be clear why the number of local degrees of freedom 
of interacting massive gravity theories does not change discontinuously in most of the 
limits in which the linearized theory is discontinuous. We should expect discontinuities 
in decoupling limits: the number of local degrees of freedom is reduced in the TMG 
limit of GMG, and the EC limit of TMG, because we are taking a limit in which one 
degree of freedom becomes inaccessible; this is not a physical discontinuity. The only real 
discontinuity occurs when the non-linear theory acquires an enhanced gauge invariance, 
and this occurs only in limits that yield 3D conformal gravity. These considerations 
explain the results of the above table. 

Another result of this paper is an alternative Chern-Simons-like form of the action 
for TMG and GMG in which the Lagrange multipliers for the secondary constraints are 
promoted to new one- form fields. In the TMG case, this action is just the sum of the 
CS actions for 3D GR and 3D conformal gravity. As for our Hamiltonian form of NMG 
and GMG, there is only a "branch-equivalence" to the original CS-like actions, but this 
just means that we have a slightly new 3D massive gravity model that is known to be 
unitary at least when linearized about one branch of its solutions. One might think 
that a time/space decomposition of this alternative action would lead to a Hamiltonian 
formulation that preserves the manifest 3-space covariance in the same way as CS models. 
However, since the space components of the new one-form fields are zero as a consequence 
of the field equations we now need new non-covariant secondary constraints, leading to 
a Hamiltonian formulation that is now "non-minimal" but otherwise equivalent to the 
minimal formulation. 

Following the work of [4j on 3D massive gravity models, progress has been made 
towards the construction of a non-linear and ghost-free 4D theory of massive gravity, 
e.g. [3 1J and [32] • The model presented in the latter paper has received much attention 
although it was, for a while, a challenging technical problem to prove that it is ghost-free; 
see e.g. [33] and references therein. In recent work [34J, this 4D massive gravity model has 
been reformulated in a vielbein language that is reminiscent of our CS-like formulation 
of 3D massive gravity models, so it would be interesting to see whether the methods that 
we have used here could also be used to simplify its Hamiltonian formulation. 
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Note added. Following the posting of the original version of this paper to the arXiv, 
a revised version of ilarXiv: 1208.0339 1 was posted, with a revised title, in which a Hamil- 
tonian analysis of the model referred to here as "massless NMG" is presented (35], with 
conclusions that are in accord with ours. 
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